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Abstract 

We discuss differential- versus integral-equation based methods describing 
out-of thermal equilibrium systems and emphasize the importance of a well 
defined reduction to statistical observables. Applying the projection operator 
approach, we investigate on the time evolution of expectation values of linear 
and quadratic polynomials in position and momentum for a statistical anhar- 
monic oscillator with quartic potential. Based on the exact integro-differential 
equations of motion, we study the first and naive second order approximation 
which breaks down at secular time-scales. A method is proposed to improve 
the expansion by a non-perturbative resummation of all quadratic operator 
correlators consistent with energy conservation for all times. Motion cannot 
be described by an effective Hamiltonian local in time reflecting non-unitarity 
of the dissipative entropy generating evolution. We numerically integrate the 
consistently improved equations of motion for large times. We relate entropy 
to the uncertainty product, both being expressible in terms of the observables 
under consideration. 
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1 Introduction 



Non-equilibrium aspects of quantum field systems are of actual interest. Among other prob- 
lems, one challenge is the description of time evolution of macroscopic quantities such as the 
expectation value of the field strength in systems out of thermal equilibrium. The problem, 
essentially the basic question for macroscopic dynamics of statistical systems, and thus inter- 
esting in itself, also plays a prominent role in context with cosmological inflationary phases. 
Although a plethora of equivalent descriptions of the collisionless limit are available, it is 
difficult to elaborate methods which allow to go beyond this thermodynamically trivial limit. 
One purpose of this paper is to promote a strategy based on which time evolution can be 
described in a consistent and systematic way. 

Although the basic concepts of thermodynamic statistical theory are known for a long 
time, there still appears to be some confusion about terminology and underlying concepts 
of non-equilibrium systems and their dynamical description. Another task of this work is to 
introduce, review and discuss some of the basic ideas in a brief but systematic way. Tech- 
nically, we adopt the projection operator method as one consistent method to describe time 
evolution of physical quantities in statistical systems. 

To keep complexity minimal, we consider the toy model of the quantum Duffing oscillator 
for practical calculations, which nevertheless is a system of sufficient complexity to exhibit 
important features of non-equilibrium thermodynamics. For systematic reasons, we discuss 
some properties of that system in the statistically trivial, entropy conserving, approximation 
first. We present criteria for permissible initial data, study static solutions and their stability, 
construct the first order integrals of motion and the first order effective Hamiltonian and 
Lagrangian. The system exhibits phenomena of parametric resonance which best may be 
displayed in the small coupling limit. 

Finally, we attack our main task and investigate on proper non-equilibrium features. In a 
first step, that requires the study of second order effects which, however, turn out to violate 
energy conservation at long time scales in the strict power-series expansion. The main result of 
this investigation is to improve the second order result by taking into account the non-unitary 
effective evolution of observables which renders the approximation scheme self-consistent and 
solves the problem of non-conservation of constants of motion. In the appendix, we explicitly 
derive the basic uncertainty-entropy relation for the set of quadratic observables. 



2 Definitions and general properties 

(i) A mixed state of a quantum system (configuration), both in zero dimensional quantum 
mechanics as well as in quantum field theory, is described by a density operator, which, in 
order to allow a probability interpretation, must be a hermitian trace-class operator with 
positive eigenvalues. It is an intrinsic feature of quantum systems that the density matrix is 
Active in the sense that only its diagonal elements correspond to physical probabilities, the 
other dependencies being phases which enter in expectation values via interference effects^. 
Alternatively, one may characterize a configuration by the expectation values of hermitian 
operators. A generic set of those representing a complete set of observables is given by the 
mutually orthogonal hermitian projectors constructed from the eigenvectors of the density 
matrix. A complete set of (not necessarily commuting) observables contains all information 

2 That is the fundamental difference to classical phase space averages. 
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about the density matrix such that any observable can be expressed in terms of the complete 
set. 

(ii) The system may be assigned a dynamical structure. Motion is defined as a sequence of 
possible states having a constant expectation value of the Hamiltonian operator, the energy 
of the configuration. Quantum mechanical time parameterizes those configurations which 
are assumed to have time-independent probabilities and phases for autonomous systems. 
The time evolution generated by the Hamiltonian can be expressed by first order differential 
operator-equations in time for the density matrix (von Neumann equation) in the Schrodinger 
representation. Hermitian Hamilton operators generate unitary time evolution which in turn 
is necessary for compatibility of the probability interpretation with time evolution. Once 
the problem of time evolution of the mixed state is solved for a given initial configuration, 
observables can be calculated as expectation values from the evolved density matrix. 

(iii) The statistical description of a system is based on a reduction procedure selecting an 
in general much smaller number of (macroscopic) observed quantities out of the complete set 
of observables. This subset defines the level of description. The process of reduction from 
the density operator to the set of the so-called relevant quantities in general involves loss of 
information. Here, we will concentrate on levels of descriptions that once chosen, will be kept 
fixed during time evolution. That constraint can be relaxed too if necessary Q. 

(iv) The reduction can be dynamically trivial if it commutes with time evolution. In that 
case of dynamical closure, the Liouvillian maps the operators of the level of description on 
a linear combination of them. The observables of a level of description at a certain time 
are sufficient to determine them at any time and evolution induces neither information gain 
nor information loss at the level of description, the associated entropy being constant. On 
account of the canonical commutation relation, the most general Hamiltonian admitting a 
finite dimensional dynamically closed level of observation can contain constant, linear and 
quadratic expressions in position and momentum operators^]. The corresponding dynamically 
closed sets of observables correspond to sums of polynomials of finite order in the canonical 
operators. If the reduction is non-trivial, one can extend the level of observation to render it 
trivial. That may involve an infinite number of operators in which case the system is a truly 
interacting one, and the only dynamically closed set of operators corresponds to the complete 
set of observables. 

(v) For truly interacting dynamical systems, the complete initial density matrix influences 
observables at later times. Its definition calls for an additional principle to construct it from 
the reduced set of initial data. Information theory proposes to apply Shannon's theory of 
entropy ||] to that physical problem. Jaynes' principle of maximum entropy [|| fixes the 
generalized canonical density operator as initial condition. It can be shown not to contain 
more information than the initial set of observables. We want to point out that this choice is 
the statistically most probable, but the underlying concept of ensemble averages of identical 
systems evolving from variant initial preparations does not imply the actual preparation of 
the system in that state. 

We confront two strategies to arrive at a description of the motion of observables. 

(A) Solve the complete problem of time evolution for the mixed initial state and extract the 
interesting quantities. The reduction procedure explicitly is postponed after time evolution. 

(B) Construct equations of motion for the relevant quantities. They turn out to be integro- 

3 Regardless of additional dependencies on other c-number quantities, including time and classical field 
strength, we call it free (effective) Hamiltonian. 
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differential equation in time necessary to include the effects of the past history on the relevant 
observables. 

Solution (A) appears to be be more attractive for field theorists since methods of field 
theory can be applied with rather little modification. The problem boils down to solve the 
initial value problem for mixed states for which the appropriate tool of the closed-time-path 
method has been developed ||]. We sketch the systematic procedure (A) according to the 
program outlined in (i)-(v), and discuss associated technical and conceptual problems. 

Starting with the equations of motion in differential form, one gets an infinite hierar- 
chy of first order differential equations for the complete set of observables, corresponding to 
the Schwinger-Dyson hierarchy of equal time correlators, and analogous to the Bogoliubov- 
Born-Green-Kirkwood-Yvon hierarchy of many-particle physics ||. The hierarchy may be 
represented in compact form in terms of a classical effective potential generating the equal- 
time correlators || by transforming to the Wigner representation. In principle, the reduction 
amounts to eliminate the irrelevant quantities resulting in a finite set of infinite order differ- 
ential equations for the reduced set of observables the integration of which involves an infinite 
number of initial conditions. At that point the statistical assumption about the initial den- 
sity matrix enters to determine all derivatives of the observables at the initial time to- The 
essential feature of those effective equations of motion is non-locality in time which becomes 
apparent when converting them into integral equations. Suppose we only knew the relevant 
observables at a later time t\ > to, then neither retrodictions nor predictions can be made 
on account of the lacking knowledge on the derivatives at t\, or the complete history of ob- 
servables. In that sense, reduction necessarily goes along with current temporal information 
loss. A complete calculation of the present observables can only be done on grounds of the 
complete past initial density matrix where missing information was induced by the princi- 
ple of maximum uncertainty. The corresponding reduced entropy at later times is always 
larger than the entropy of the initial state which introduces thermodynamic irreversibility in 
a natural way. 

We want argue that, strictly speaking, the cut hierarchy is not a permissible dynamical 
approximation. There, one keeps only those equations which involve the time derivatives 
of the variables of interest. In order to arrive at a closed set of equations, the irrelevant 
quantities generally still being present in the system have to be expressed somehow in terms 
of the relevant ones. There is no systematics of how to do that. In some particular cases, 
as the large N-expansion of O(N) symmetric scalar theories, the hierarchy by chance closes 
without further assumptions |^]. In the classical theory of gases Boltzmann introduced the 
famous Stofizahlansatz for two particle correlators at that point. Generally dynamics then is 
represented by a finite closed set of effective (macroscopic) equations of motion. Of course no 
statistical assumptions about the complete initial state has to be inferred. 

We want to point out that the truncation of the differential hierarchy explicitly alters 
the dynamics of the system in an uncontrolled way. Any finite order hierarchy of equations 
of motion exhibits features of a closed dynamic system discussed in (iv) and is thus a pri- 
ori inequivalent to the reduced truly interacting evolution. It is always local in time which 
presumes memory loss at microscopic scales and automatically reduces to a Markovian de- 
scription, even in cases where the complete system behaves differently. Some truncations are 
compatible with the existence of a macroscopic effective Hamiltonian (Lagrangian) being a 
function of a finite set of effective variables. The very existence of that Hamiltonian descrip- 
tion, however, automatically conserves the relevant reduced von Neumann entropy such that 
the system can neither exhibit equilibration nor thermodynamic irreversibility in the strict 
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sensed). Also, due to the Liouville-theorem the volume of cells in the phase space remains 
constant under effective Hamiltonian evolution. The volume, on the other hand, corresponds 
to the classical particle number density, which enters in the entropy functional of Boltzmann. 

The argument may best be illustrated for the class of effectively Gaussian approximations 
which can be obtained by time dependent Hartree-Fock variational methods Q , the method 
of equal-time Green functions |^, flCfl ) optimized expansions [11] or the truncation of the 
Schwinger-Dyson hierarchy to the one and two particle Green function. All those methods 
are physically equivalent in the Gaussian limit where the Hamiltonian is assumed to have free 
form, or, equivalently the density matrix is approximated by a generalized Gaussian wave 
package. One considers time dependent spatially bilocal functions — essentially equal-time 
two point correlators — as coefficients in the quadratic expressions which are determined by 
the corresponding variational and expansion methods. But those bilocal coefficients enter in 
the generalized uncertainty which is a constant of time in the Gaussian approximation. On 
account of the uncertainty-entropy relation, also entropy remains constant contradicting in- 
formation loss inherent in the incomplete description. Under certain conditions, the effective 
Gaussian theory can even be shown to have a representation in terms of a free theory after 
requantisation [12|. Although this approximation may supply a fairly good description at 
short terms or close to equilibrium, they are thermodynamically trivial approximations and 
cannot account for full long term evolution since they effectively correspond to the collision- 
less^] limit of the physical system. We will show that in fact an effective Gaussian density 
matrix, whatever its time dependence be, does not evolve to the Gaussian approximation of 
the interacting density matrix at large times. The 'equilibration' observed in those approx- 
imations [l3|] can thus not be identified with physical entropy generating processes, but 
converges to static values by dephasing effects. 

In general, the dynamics of a cut hierarchy effectively corresponds to current preparation 
of the system to have the canonical density matrix of the level of observation. Thus, dynamical 
closure is attained at the price of effectively time dependent probabilities being typical for 
non-autonomous systems. The procedure presumes that the density matrix decoheres for 
macroscopic observables including the existence of a classical limit for the variables of the 
level of observation. Such an effective theory may exist, depending on the particular system 
under consideration, but it remains to show that it represents motion also at macroscopically 
large time scales. 

The problems with the differential representations can be avoided if we adopt solution 
(B) which intrinsically accounts for non-locality in time. The practical method we apply 
is the projection operator method in Schrodinger representation which yields the Robertson 
equation |Q being equivalent to the generalized Langevin-equation ||, |l5| in the Heisenberg 
representation. Although those equations of motion are closed in the variables of the level of 
observation only, expectation values of arbitrary operators can be expressed by non-local time 
dependence. The exact reduced equations of motion turn out to have a rather complicated 
structure and suitable approximation schemes have to be developed within that framework. 
Splitting off a dynamically closed part from the Hamiltonian, the entropy generating contri- 
butions involve quadratic and higher powers in the remaining truly interacting Hamiltonian. 

4 Effective irreversibility, sometimes misleadingly related to ergodicity Q , may occur in infinite-dimensional 
systems ||, but has to be justified a posteriori and cannot be concluded on grounds of an approximation 
scheme. 

5 Strictly speaking, the term collisionless is not precise since the limit does account for effective elastic 
scattering. 



5 



There, a non-unitary operator enters, responsible for time evolution in the projected sub- 
space. Approximating that operator by a unitary one again amounts to introduce a Gaussian 
approximation which differs formally only at third order from the exact solution but never- 
theless involves unphysical consequences. It conflicts with entropy generated at second order 
and, more severely, results in non-conservation of the energy at secular scales. 

3 Basic Setup 

We briefly sketch the result of the projection operator approach [l5|]. For convenience, 
we define a level of observation by the finite set of hermitian operators £ = {J~v} including 
To = 1 . For the corresponding operator expectation values g u (t) = tv (Tu p{t)) , the closed 
exact equation of motion is found to read 

^-tr (.F„ft(*)) = -i tr {T V L o K(t)) - f dt' tr (T V L o T(t, t') o Q(t') o L o K(t')) . (1) 
at Jo 

The accompanying canonical density operator 1Z(t) = exp(— Hv(f)J~v) which minimizes the 
entropy in £ contains time dependent Lagrange multipliers /j,(t) which are functions of g v (t) 
such that g u {t) = tr (JFjTZit)). We emphasize that TZ(t) does not evolve with L and must 
not be confused with the density-operator p(t) of the system. All trace expressions at the 
r.h.s. of (||) can at least in principle be expressed in terms of the g v , and the system is a 
closed integro-differential equation in the c-number expectation values. The projector Q(t) is 
defined by 

tr(OxQ(t) o Q 2 ) = tx{0 1 2 ) - ^f 1 ^ tr(F u 2 ), (2) 

and the non-unitary evolution operator T(t, t') is a solution of 
d d 

—T(t,t') = -iQ(t)oLoT(t,t') and — T(i, t) = iT(t, t') o Q(t) o L (3) 

with initial condition T(t,t) = 1. 

Here, we will complete £ such that the action of the Hamiltonian [H, X] = L o X can 
be split into L = Lq + L\, and the free part be dynamically closed with respect to the level 
of observation, Lq o T v = [Hq,^] = ^l^^J 7 ^. In that case, in the integral (|l|), the complete 
Liouvillian can be replaced by L\. 

Expectation values of operators which are not in the linear hull of £ get additional con- 
tributions to their 7£-averages, 

tr(0p(t)) = tr (OR(t)) - i [ dt'ix (OT(t, t') o Q(t') o L o K(t')) (4) 

J o 

where again the dynamically closed part in L does not contribute to the integral. The integral 
vanishes for O G £. 
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4 Zero-dimensional system 



We will investigate on the time evolution of an anharmonic oscillator defined by the Hamil- 
tonian 

H = H + H 1 , H = ^(p 2 + sx 2 ), H 1 = ^x 4 (5) 

with canonical commutator [x,p] = i. The parameter s = ±1 allows to switch to the broken 
symmetry case. We define the level of description by a dynamically closed set of observables 
for Hq given by £ = {1 , x,p, x 2 ,p 2 , xp + px}. That set plays a preferred role since, due to 
its bilinear nature, the accompanying density operator lZ(t) is quasi-Gaussian which allows 
to apply a modified Wick theorem in the evaluation of expectation values. Inclusion of the 
quadratic variables also ensures the existence of the accompanying canonical operator, and 
allows for direct comparison with approximations of Gaussian type. The equations of motion 
are found to read 

f t x = p (6) 
^-p = -sx + 4Ax 3 - 6Xxx 2 + S p (i) (7) 

Jt X2 = W (8) 
d __o 

— p2 = 8Xpx — sw — 6Xwx2 + S p 2 (t) (9) 

j t w = 2p 2 + S\x A -2sx 2 -l2\x 2 2 + T lw (t) (10) 

with x = (x) , p = (p),X2 = (x 2 ), P2 = (p 2 ) , w = (xp + px) and 

= - /* di!tv(T v L x o T(t, t') o Q(t') oL l0 K{t')) (11) 
Jo 

accounting for entropy increase. Note that x and ±2 do not get further corrections since the 
potential is a function of position only. The energy is of the form 

tv(Hp(t)) = i(p 2 + sx 2 ) + ^(3xi-2x 4 )+e(t), 

e(t) = -i f l dt'ti(H 1 T(t,t')oQ(t')oL 1 oTl(t')), (12) 
Jo 

and is conserved in time. The key challenge is to find a sensible approximation for the operator 
T(t,t r ). However it is worthwhile to study the first order system without integrals first. 



5 First order approximation 

5.1 Permissible initial conditions 

The first order system may be integrated numerically with initial conditions for x, . . . , w. 
Those as well as their time evolved counterparts are subject to the positivity conditions, 
%2 > 0,P2 > 0, 5 2 x = X2 — x 2 > 0, S 2 p = P2 — P 2 > 0, and the uncertainty relation 
(S 2 w = w — 2xp) 

A 2 :={5 2 x){5 2 p)-\{5 2 W ) 2 >\. (13) 
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which follows immediately from Schwarz' inequality relation^ (A 2 ) (B 2 ) > \ (AB) \ 2 . Coeffi- 
cients that do not fulfill that relations cannot appear as expectation values averaged with a 
physical accompanying density operator 1Z. 

5.2 Static solutions and their stability 

Static solutions of (||-[H]) are found to read 



x = ±y + 4 6 A AX2 , p 2 = -^{s + 2\x 2 ){s + 6\x 2 ), w=p = 0. (14) 

In the unbroken case s = 1, for any choice of x 2 the value p 2 is always negative. We conclude 
that whatever motion be, for permissible initial conditions the system never attains a state 
where all observables become time independent. 

In the broken case s = —1, however, the positivity requirement for p 2 permits static 
solutions in the range 1/(6A) < x 2 < 1/(2 A). They may have expectation values for position 
with |x| < l/\/2Aj where the bounds correspond to the classical local minima of the potential. 
Classically, one would not expect static solutions for values of x other than the minimum 
l/\/2A. An additional restriction is given by the uncertainty relation, 

(1 - 2Ax 2 ) 2 (-l + 6Ax 2 ) > 4A 2 . (15) 

For A -C 1, this condition narrows the region of static solutions 1/(6A) + 3A/2 + 0(X 2 ) < 
x 2 < 1/(2A) — l/v2 + by a small amount. For increasing A, the window gets smaller 

and finally completely vanishes at A = §y§ ~ 0.18, corresponding to x 2 = \\J\ ~ 1-53 and 

x = (27/2) 1,/4 /2 ~ 0.96. For larger couplings, no static solutions exist. In that case, the 
minimal quantum uncertainty suffices to overcome the potential barrier and the motion can 
escape the initial half. 

To study the stability of the static points in the broken phase, we consider a set of 
initial conditions with p = w = and x 2 ,p 2 satisfying the staticity requirement. We chose 
x(0) = x s ^ a ^ lc +<5, where 5/x <C 1 to model perturbed initial conditions. For 5 > a numerical 
integration shows that the position observable periodically increases to a maximum, at which 
the width 5 2 x develops a local minimum. The motion of x remains in the initial half, but the 
amplitude is large with respect to the initial coordinate. 

For 5 < 0, the qualitative picture is completely different. There the average position 
performs an oscillation to the negative initial value, while the width oscillates only slightly. It 
is surprising that motion is just opposite to the classically expected roll-down in the potential. 
Again quantum effects largely dominate over the classical picture. 

In both cases, the motion of the position coordinate does not remain near the initial 
preparation. In that sense, static solutions are not attractive, and cannot be regarded as 
stable. 

5.3 First order integrals of motion 



Rewriting the system (|6-1C) it in terms of the quantum widths X = \'o 2 x, P = ^/fPp, W 2 = 

^W 2 = 2P 2 -2(s + 6\x 2 )X 2 , ^P 2 = -( s + 6\ X2 )W 2 , ^X 2 = W 2 , (16) 

dt dt dt 

6 This condition generally gives stronger constraints than the usual procedure to minimize the positive- 
definite expression ((A — AB) 2 ) > 0. 
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one verifies the time independence of A. It is remarkable that the uncertainty appears as 
constant of motion which moreover in independent of the coupling. Physically, on account of 
the uncertainty-entropy relation (see Appendix), constant uncertainty means that quantum 
effects do not increase entropy. To first order, the value of the set of observables £ at a 
given time is sufficient to determine their time evolution locally in time, and to express all 
observables of the system. 

A second constant of motion is given by 

h = \P2 + \sx 2 + ^(3x 2 - 2x 4 ) (17) 

which equals the expectation value of the Hamiltonian averaged with 1Z. That statement is 
not so trivial as seems at first sight. It is a consequence of the fact of dynamical closure of 
Hq which makes the corrections in the integral in (|l2|) of second order. 

Let us further make a comparison with the standard perturbative approximation. There, 
one expands H around the classical expectation value of momentum and position, i.e. x = 
x + 5x, where 5x is supposed to be small. Applying that approximation to the Hamiltonian 
and keeping only quadratic terms in Sx, one finds for the A-dependent contribution to the 
energy X(Qx 2 x 2 — 5x 4 )/2, which differs from h if 5x is not small, to wit, if quantum or 
statistical corrections become important which is in fact the case. Note also that in the case 
of symmetric configurations x = 0, the quartic potential part does not contribute to energy 
in that approximation at all. Thus the standard perturbative first order effective potential 
does not coincide with the consistent Gaussian dynamical potential, in agreement with [||. 

5.4 First Order Effective Hamiltonian and effective Lagrangian 

The energy integral can be used to eliminate p,p2,w from the system (|6 10 ) , 

x = — sx + 4Ax 3 — 6Xxx 2 , (18) 
X2 = 4h - 4sx 2 - I8\xl + l2Xx 4 . (19) 

This system is not a Hamiltonian system in the kinetic momenta x, x 2 . A necessary condition 
of the existence of a Hamiltonian x 2 + x\ + J(x, x 2 ) would require J to be a solution of 

— = -x(x,x 2 ), — = -X2(x,x 2 ). (20) 

However, the integrability condition d 2 J/(dxdx 2 ) = d 2 J/(dx 2 dx) is not fulfilled by the cor- 
responding expressions of (18], 19). 



Nevertheless, an effective Hamiltonian which reproduces the equations motion of (|-10) 
together with the constant of uncertainty A 2 can be constructed. We introduce the canonical 



pairs (x,ir) and (x 2 ,tt 2 ) and cast the conserved energy (|17| ) in canonical variables. Then, p 2 
becomes a function of (x, x 2 , n, tt 2 ) and has to be consistent with the Hamiltonian equations 
of motions, 

dH 1 dp 2 . dH 1 dp 2 

P = X = ^ = 2^ w = X2 = d7 2 = 2d7 2 (21) 

On the other hand, the constant uncertainty requires p 2 to be a solution of the partial differ- 
ential equation 

i (d P2 \ 2 -^P2dp 2 + x 2 fdp 2 \ 2 =p2x2 _ A2 (22) 



16 \d7r 2 J 4 dir 2 dir 4 \ dir 
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Splitting off a summand A 2 / X 2 the remaining function has to be a homogeneous expression 
in 7r, 7T2 of order two. We find 



A 2 

P2 = + (4X 2 7T2 +4x7r 2 7r + 7r 2 ), p= (2x7T 2 + 7t), (23) 

which relates observables to the canonical variables. The corresponding Hamiltonian function 
has the form 

H(x, 7r; x 2 , vr 2 ) = 2x 2 7r 2 + 2x7r 2 7r + ^7r 2 + V(x, x 2 ), 

2F(x,x 2 ) = — + sx 2 + A(3xi - 2S 4 ) (24) 

It can be checked that the remaining Hamiltonian equations of motion dH/dx = —tt, dH/dx2 = 
— 7r 2 are in fact compatible with the original set of equations of motion. The corresponding 
effective Lagrangian is related to H by a Legendre transformation, L = 7r 2 i 2 + ttx — H . Eq. 
( |2l| ) supplies the necessary relations of canonical and kinematical momenta, and we finally 
get 

L = ^X 2 + h 2 -V(x,x 2 ). (25) 

In terms of the coordinates x, X, the Lagrangian is composed of a kinematical kinetic contri- 
bution, a potential (si 2 + sx 2 + A(3X 4 + 6X 2 x 2 + x 4 ))/2 and an additional term _4 2 /(2X 2 ) 
having the form of the energy of circular motion with angular momentum A. 

5.5 The small coupling limit and parametric resonance. 



In the strict limit A = the equations of motion ([L^jl^) decouple trivially. However, that limit 
does not provide a sensible starting point for a perturbative expansion since the perturbation 
series turns out to be non-analytic around zero coupling. We thus concentrate on the limit 
A — » + in which the equation of motion for x has the structure 

x- + xuj 2 {t) = (26) 

of a linear oscillator with a time dependent frequency- factor. It can be seen by completely 
eliminating x 2 from the system (18,1^) that for finite x the quantities x 2 ,p 2 and w consistently 



have to be taken of order 1/A. Consequently, also the energy being a function of those 
quadratic expectation values grows with the inverse coupling and will be rescaled into the 
finite quantity h = 18Xh + 1. The time dependent mass-term uj 2 = (s + 2z\//i)/3 can be 
rewritten in terms of the function z = (9Ax 2 + s)/Vn which is a solution of z" + z 2 — 1 = 0, 
where the prime denotes differentiation with respect to the rescaled time r = t(2\/7i) 1 ^ 2 . 
Note that in this limit, there is no back-reaction on the frequency-factor by the motion of the 
position coordinate. A first integral is given by 

l -{z') 2 =A + z- Z ^ (27) 

where A encodes the initial conditions of x 2 . Solutions for z(t) are periodic if \A\ < 2/3 with 
minima in the strip \z\ < 1 and maxima for 1 < z < 2. 
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Figure 1: Time evolution of the effective mass 

uJ 2 (t) (dashed line) and x(t) with parameters ^ nie 2: The same as in Fi §- but with 
h = 3.1, z = -0.56, x = 1, x' = z' = 0. h= 1M ■ 



For physical solutions of the system (£6-27) we have to keep in mind that only positive 
values of X2 can be generated by hermitian density matrices^. Consequently, zvh — s has to 
remain positive in the course of time evolution. In the case of the unbroken potential s = 1, 



this further restricts motion of z to the strip 2 < 2zVh < (— 1 + y 6 + 3h) where the bounds 
correspond to the local extrema of a periodic motion during which u 2 remains positive. 

The more interesting case is the broken phase s = — 1. There, the solutions which fulfill 



the positivity requirement for X2 can now be found in the domain —2 < 2zVh < (1 + y 6 + 3h) 
and correspond all to oscillatory z. However, if the initial conditions are varied, the behavior 

of uj 2 becomes quantitatively different. In particular, in the range 1 < 2z\fh < (1 + y^6 + 3h), 

the factor uj 2 is positive but becomes negative for —2 < 2z\fh < 1. In the first region, the 
motion of x is quasi-harmonic with variable frequency, but turns into exponential behavior 
if the lower strip of z is reached by time evolution. Physically, oscillatory and tunneling 
phases alternate. Negative square of frequency flips the sign of the curvature of x and bends 
trajectories in a direction opposite to the oscillations. However, since there is no back reaction 
on the motion of x from the frequency- factor, the tunneling phases are short, and exponential 
growth stops before the numerical value of x can become large. The period of x is roughly 
determined by to 2 averaged over one period of z. Depending on the quotient of those two 
periods, a variety of resonance phenomena can appear. Figure (1) displays the case of a 
frequency coefficient close to 1 : 2. The periodic kick of the motion of x results in large 
resonant amplitudes. Figure (2) shows a frequency quotient of 3 : 8. The dents in the motion 
of x correspond to tunneling phases. 



7 8 2 x > is trivially valid since Xx2 is of the same order as x. 
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6 Second order Motion 



Realistic thermodynamic behavior includes entropy generating processes which necessitates 
the inclusion of the non-local integral contributions in (|l|). The evaluation of that term, 
however, requires an approximation for the evolution operator T(t,t'). A strict expansion 
in powers of the coupling constant to second order amounts to the replacement L — > Lq in 



the definition of T(t,t') which has been applied recently in the context of field theory [17]. 
On account of the dynamical closure, T(t, t') even exponentiates to the product of the free 
unitary time evolution operator with Q. That approximation, however, suffers from non- 
conservation of the energy ([12; ) at secular scales. A numerical integration shows that the 
runaway of energy sets in even at times before the interesting transition from regular motion 
to the stochastic entropy generating phase appears. Moreover, since energy turns out to 
increase even exponentially, the expectation values also grow in an unphysical manner. A 
first attempt to cure the situation is to replace the free Hamiltonian in T(t, t') by a resummed 
quadratic one. At first sight this concept appears promising, since it still benefits from a 
dynamic closure relation. We discuss this approach first before we turn to a further essential 
modification. 

6.1 Resummed Hamiltonian approach 

The resummed first order Hamiltonian is defined to be of at most quadratic in position 
and momenta but with explicit time dependence through coefficients (a, ... ,e), and has to 
reproduce the equations of motion (P-^0|). We make the Ansatz H r (t) = Hq + AiJ[(i) with 

H[{t) = ap + bx + c(xp + px) + dx 2 + ep 2 . (28) 

Plugging this expression into ([!]) and keeping terms to first order in the coupling, the equations 
of motion turn into 

k = p(l + 2Ae) + 2Xcx + Xa, p = -x(s + 2Xd) - 2Xcp - Xb, 

x 2 = w(l + 2Xe) + 2Xax + kXcx 2 , p2 = ~w(s + 2Xd) - AXcp 2 - 2Xbp, (29) 

w = 2p 2 (l + 2Ae) -2x 2 {s + 2Xd) -2Xbx + 2Xap, 

which when compared with the original first order set of equations defines 

a = t,(2px 2 - wx), b = -4x 3 + £(2p 2 x-pw), c=^W 2 , d = 3x 2 -£P 2 , e = -£X 2 . 

(30) 

The one-parameter solution parameterized by £ is a consequence of time independent un- 
certainty which still is an integral of motion of Eq. (|29|). The remaining freedom has to be 
fixed by the condition that H r is a constant of time which is equivalent to the fact that the 
expectation value of W equals the expectation value of H, 

c = 3(2x*-x 2 ) ^ 

W 2 — AX2P2 + Q{X2P 2 + X 2 P2 — WXp) 

Resummation is achieved by absorbing the effective Hamiltonian into the new free one, 
and compensating by an interaction counter term, 

H = H r + H 2 , H 2 = H 1 -XH r 1 {t). (32) 
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In the integrals of the dissipative contributions in ([ll]), we have to replace L\ by L2. However, 
the set of observables £ is still dynamically closed with respect to H r such that only Hi of 
H2 contributes. Finally, the resummation enters in the expansion of the evolution operator 
T(t,t') where we replace L — * L r (t) in (^). We have to consider the particular combination 
T(t,t') o Q(t') appearing in the entropy generating contribution which can be integrated 
formally directly from its definition, 

T r {t,t') o Q(t') = Texp (-iJ drQ(r) o L r (r) o \ Q(t'). (33) 

At that point, dynamical closure guarantees the important relation Q(t") o L r (t) o Q(t') = 
L r (t) o Q(t') valid even in the case where the time arguments of the projectors differ. All Q 
but the one at the most r.h.s in the summands of the formal exponential can be dropped, and 
we arrive at the desired relation 

T r (t,t')oQ(t') =U(t,t')oQ(t'), U(t,t') = Texp(-iJ*dTL r (T)o\. (34) 

We want to point out that this relation does not automatically imply the equality of T r (t, t') 
and U(t, t'). The replacement is non-trivial also in the sense that acting on the projector, the 
action of a non-unitary operator is replaced a unitary one. 

It turns out to be useful to let IA act on the operators to its left, tr ((lft(t, t') o 0\) o 2 ) = 
tr (0i oU(t,l/) o O2). On account of the two properties of L r being a derivation and anti- 
hermitian, the adjoint action simply amounts to exchange time arguments, W(t, t') = Hit' , t). 
We introduce the evolved operator 0(t,t') = W(t,t') o O which solves 

d t ,Q(t, t') = -iL r {t') o 0(t, t') (35) 

and is related to the Dirac representation induced from U by Oz>(t, t') = U(t, t')oO = 0(t', t). 
Since the Liouvillian is a derivation, it suffices to investigate on the evolution of x,p, l 
representing an operator basis. Dynamical closure guarantees the success of the Ansatz 

x(t,t') = xa x (t,t') +p(3 x (t,t') +-/ x {t,t'), p(t,t') = xa p (t,t') +p(3 p (t,t') +j p (t,t') (36) 

with a x (t,t'), . . . , 7p(t, t') to be determined by coefficient comparison. That gives rise to the 
system 

d t ,a x (t,t') = -2Xc(t')a x {t,t') + (s + 2Xd(t'))p x (t,t') 
d t ,p x {t,t') = -(l + 2\e(t'))a x (t,t') + 2\c(t')(3 x (t,t') 
d t Hx(t,t') = -\a(t')a x (t,t') + \b(t f )p x (t,t>). (37) 

Integrals are found to read 

a x {t,t')x(t') +(3 x (t,t')p(t') +7*(M') = x(t), (38) 

which corresponds to the quantity tr(x(t, t')H(t')) independent of t' , and 

a 2 x (t, t')X 2 (t') + 0l(t, t')P 2 (t') + a x (t, t')(3 x {t, t')W 2 (t') = X 2 (t), (39) 
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expressing the if independence of tr(x(t,t')x(t,t')TZ(t')). We remark that these constants 
imply unitary evolution of the accompanying operator, U(t,t')lZ(t') = lZ(t), and are com- 
patible with a time independent uncertainty product. The complete solution of ( p7| ) can be 
parameterized by the angle variable 

y(t, t') = J* drA (x- 2 {t) - 2Af (r)) . (40) 

We find 

a x (t, t') = A~ l X(t)P(t') coe(tf (t, t') + T]{t')), (3 x {t, t') = A~ l X(t)X{t') sin(¥(t, t')) (41) 
with 

cosr](t) =AX- 1 (t)p- 1 (t), smr](t) = W 2 (t)X- 1 (t)p- 1 (t)/2, (42) 

which satisfies the initial conditions a x (t,t) = 1, /3 x (t,t) = 0. The coefficients for p(t,t') are 
also solutions of the differential equations fl37|), but with integration constants p(t) and P 2 (t) 
at the r.h.s. of (38,j39). The corresponding angle variable gets shifted to ty(t,t') — rj(t) + ir/2 
and we get 

a p (M') = -A^P^Pit^sm^it^^-^+f]^)), 

P P {t,t') = >l- 1 P(t)JSf(t / )cos(*(i,i / ) -r?(i)). (43) 



Now all necessary tools are at hand to evaluate the integrals in (11) and (12) necessary 
solve the integro-differential equation of motion in a resummed Hamiltonian approximation. 
We have integrated the system numerically and still find non-conservation of energy at typical 
time scales of order 2ir/ (4A) for A< 1. The factor 2ir corresponds to the natural period of the 
oscillator with unit circular frequency, and 1/4 can be identified to coming from the fourth 
powers of the harmonic functions in a, (5 due to the quartic interaction Hamiltonian. At a 
technical level, the problem is that quadratic Hamiltonians, even with time dependent coeffi- 
cients, automatically conserve the uncertainty A 2 , and motion stays within this equivalence 
class fixed by the initial conditions. Plugging that approximation into the entropy generating 
integrals of the equations of motion, the dissipative corrections, although formally of second 
order, do not remain small at large time scales. The effective unitary evolution operator gen- 
erated by an effective quadratic Hamiltonian cannot be used as a sensible approximation in 
the entropy generating corrections. The inconsistency becomes apparent when one compares 
the uncertainty A determined by the initial conditions with the actual uncertainty expressed 
by the evolved value of the observables. 



7 Consistent effective approximation. 

Although the effective Hamiltonian constructed in the previous section turns out to be in- 
consistent with time-variant uncertainty and time evolution at large time scales, it still ap- 
proximates time evolution locally in time at short time scales. It can be seen from numerical 
integration that the quadratic Hamiltonian adjusted to a set of given quantities x(t), . . . , w(t) 
considered as initial conditions evolves them quite accurately. Local in time, motion ap- 
pears to be generated by a quadratic Hamiltonian which can be used to approximate only 
the enveloping curves of the true trajectories, but the locally best approximating effective 
Hamiltonian does not evolve to the enveloping Hamiltonian being the best approximation at 
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a later time. One may argue that a non-quadratic extension may solve the problem of en- 
ergy conservation. However, on account of the basis of observables chosen, the accompanying 
density operator still remains quadratic, and the exact Wick theorem valid, such that the 
non-quadratic extension effectively boils down to the quadratic effective theory. 

We thus abandon the query for an effective Hamiltonian completely in favor of the con- 
struction of an effective time evolution operator. Physically, we have to admit effective equa- 
tions of motions for the observables which are not generated by a Hamiltonian function. This 
step is also necessary to account for irreversibility. 

The complete effective time evolution operator V(t, t') be subject to the following assump- 
tions. 

(a) Transitivity, V(t,t') o V(t',t") = V(t,t") 

(b) Initial value, V(i, t) = 1 . 

(c) Its adjoint acts like an exponentiated derivative operator on products, i.e. V(t, t'y o 
(0102) = (y (t>t f ) ° 0l)O^(t, ° 02)- Time evolution commutes with forming operator 
products. 

We will further specify to the 'collisionless' approximation of V(t, t') by the additional 
assumptions that 

(i) — the operator linearly maps the basis x,p,1 onto itself. That replaces the notion of 
'free motion', and generalizes dynamical closure to non-Hamiltonian evolution. 

(ii) — V(t, t') is compatible with time evolution of the accompanying density operator, i.e. 
V(t,t') o lZ(t') = TZ(t). That condition actually defines V consistently with the equations of 
motion, and the evolution of 1Z, which is just contrary to the strategy of the previous chapter, 
where the evolution operator generated by the effective Hamiltonian was responsible for time 
evolution of 1Z. 

At a technical level, time dependent uncertainty A(t) is introduced into the coefficients 



a,P which are no longer subject to the set of differential equations (37). If we define the 
evolved operators by 0(t,t') = V ] {t,t') o 0, the basis still can be expanded as in (pq). The 
requirements (ii,c) automatically lead to the relations 

tr{x(t,t')K{t')) = x(t), ti{x{t,t')x{t,t')K{t')) = x 2 (t), , . 

tr(p(t,t')K(t')) = p(t), tr(p(M>(M'W)) = P2(t), 1 j 

which we identified as constants of motion of the resummed evolution equations. In addition 
to these, we find 

tr((x(t, t')p(t, t') + p(t, t')x(t, t'))K{t')) = w{t) (45) 

the r.h.s. of which was fixed by initial conditions in resummed Hamiltonian dynamics. These 
equations are sufficient to express the evolution coefficients of x and p in a manner analogously 
to (|4l],[43|) together with (f42|), but with A replaced by the actual uncertainty A(t') in (41,1 



and A(t) in (42). The angle variable *&(t, t'), however, remains undetermined for the moment. 

In the evaluation of T(i, t') o Q(t'), we consistently have to replace the resummed unitary 
operator U(t,t') by V(t,t'). The dissipative contributions T< u (t) and e(t) are now found to 
read (^ = ^(t, t'), x = x(t), x' = x(t'), . . .) 

£„(t) = A 2 f* dt'a u {t, t'), e(t) = A 2 /' dt'e{t, t'), (46) 
J o Jo 

with 

a p (t,t') = 6^ 3 X 3 A w3 sin^(l2^ ,2 cos 2 ^-sin 2 ^)s / , 
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Figure 3: Time evolution of the uncertainty Figure 4: Time evolution of the uncertainty 
of the position operator X(t). of the momentum operator P(t). 

a P2 (t,t') = 12A'- 3 X' 3 X 2 P x 

(XX' (^4^l /2 (2cos(2* - 77) - cos rf) cos 2 * - (2 cos (2* - rj) + cos rj) sin 2 <&) 

+ xx' (§A' 2 (3 cos(2* - rj) - cos rj) cos ^ - 3 cos(^ - rj) sin 2 
+ 2pa p (t,t'), 
a w (t,t') = 48^ 3 X 3 X /3 sin* x 

[XX' (aA' 2 cos 2 * - sin 2 cos * + xx' (l2^ /2 cos 2 * - sin 2 , 
e(M') = -a w (t,lf)/8. (47) 

The last equality is due to the particular form of the operator xp+px which, when commuted 
with a homogeneous operator, evaluates to twice the degree of homogeneity in x. 

A complete solution involves also to determine the unknown angle variable ^f(t,t'). Anal- 
ogously to the construction of the effective Hamiltonian, we exploit the remaining freedom to 
require the energy to be a conserved quantity within our approximation scheme. Differenti- 
ating the corresponding expressions ( |i~2| ) with respect to time, on finds that the term de(t)/dt 
has to be compensated by the integrand of the term £„2(i)/2. That gives rise to a differential 
equation for VP with solution 

9(t,lT) = f "' <1tA{t)X- 2 {t). (48) 
Jf 

It is a non-trivial result that the angle parameter can be expressed in terms of an integral 
involving the history of A and X. That nevertheless includes all contributions which can 
be resummed in an effective evolution operator satisfying conditions (i) and (ii) in a self 
consistent way. Note also that the first order correction of (4C) can be absorbed into the 



time dependent uncertainty completely. Moreover, since uncertainty is directly related to 
entropy (see Appendix), the approximation consistently accounts for entropy variations but 
by construction keeps the energy constant in time. 

The complete set of integro-differential equations can now be evaluated numerically. A 
typical time evolution is displayed in Figs. (3-6) for the initial conditions x(0) = 1, X2(0) = 
P2(0) = 2, p(0) = w(0) = with A = 0.1 and symmetric potential s = 1. 

It can be proven in the general setting that motion is such that the entropy is always larger 
or equal than its initial value. We find that our approximation respects that property, but 
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fluctuations at intermediate times are possible. There is no sign that the system approaches 
a static point at large time scales which is not to be expected in general. 

The numerical strategy to solve the consistency problem for the integro-differential equa- 
tion was the following. Suppose we knew the exact set of variables in the time range (0,t). 
Then we benefited from the enveloping Gaussian dynamics and predicted the evolution by 
integrating the first order system with the initial data at t for a time interval (t, t + St), where 
St is typically of the order of a period of x. 

That prediction can be improved iteratively for (t,t + St) by plugging it into the integral at 
the r.h.s. which represents an inhomogeneity for the differential equation. Note that the inte- 
gral expands over the complete history of the observables which poses a problem in numerical 
integration over large time scales. Luckily, the integrand can be split into a sum of products 
of functions depending solely on t or t', such that the integrals can be reused in the extension 
procedure of time range. For times including a large number of quasi-periods of the variables, 
the contribution of the integral, although evaluated exactly, increasingly exhibits stochastic 
behavior. The method proposed thus supplies a tool to exactly calculate quasi-stochastic 
dissipative behavior which can be used to test a priori assumptions about stochastic forces. 

8 Conclusion and Outlook 

We proposed a method to calculate the time evolution of observables beyond the effective 
Hamiltonian approximation in a self-consistent systematic way thus accounting for real dis- 
sipative processes. An application of this approximation scheme to field theories is planned. 
In particular, we are interested in the influence of real collisional processes on long term 
evolution, and expect to be able to decide whether dephasing or proper entropy generating 
processes dominate in infinite-dimensional systems. A generalization to gauge and fermion 
systems is planned, and we intend to develop a systematic expansion to higher orders. Ap- 
plications may also include dynamical problems in cosmology such as inflationary phases. 
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A Entropy, Uncertainty and Effective Lagrange Multipliers 



We relate the Lagrange multipliers \x v and the relevant entropy S = \i v g v to the set of 
observables x, p, X2,P2, w. One considers 

[T P ,TZ] = -n f 1 dxn x [F^F a \R}~ x , (49) 
Jo 

which has vanishing trace. That gives rise to the set of equations /j, (T g p T^' T = 0, [.F^-Tv] = 
Tp^Tp which are not independent, but lead to 

fJ L x = V (xP 2 -±pW 2 ), ii v = r,{pX 2 - \xW\ » x 2 = - V P 2 /2, f i p 2 = -r ] X 2 /2 (50) 

where rj = ^ W /W 2 . Thus, one more relation is called for to determine the remaining un- 
known. To that end, we consider the eigenvectors Q± of \oglZ being linear combinations 
of the operators T v . The eigenvalue equation [£/±,log7£] = i±G± implies ti:{Q T Q±R,) = 
e^ ± ti{Q±Qz^TV), with the conjugate pairs of eigenvectors one of which is found to have the 
eigenvalues £± = ±2^/ \i x iii p i — ^ = ±r/A. On the other hand, the traces can be evaluated 
using Wick's theorem, and one finds 




which expresses the remaining unknown rj in terms of the observables. From the constant 
normalization tvTZ(t) = we get /to + fngi = 0, and thus S = mgi. The summands can be 
combined in terms of A(t) and we finally get 

S(t) = (A(t) + i) log(,4(t) + ~) - (A(t) - I) \og(A(t) - i) (52) 

where a possible integration constant was chosen such that entropy vanishes at the minimal 
possible uncertainty corresponding to A = |. S(t) is a monotonically increasing function of 
A(t). This relation was found in a different context earlier |fl6|] . 
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